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/ Example for ratiometric approach to noise cancelling */ V 
/define LAMBDA 0.95 y ' 



void OxiLSL_NC( int 
int 
int 
int 
int 
int 
int 



reset, 

passes, 

*signal_l, 

*signal_2, 

*signal_3, 

* targe t_i, 

*target_2) 



{ 



int 

static int 
static float 



if ii, k, m, n, contraction; 
*s_a, *s_b, *s_c, 
Delta_sqr, scale, noise_ref; 



*out_a, *out_c; 



if( reset » 


TRUE) { 


S a wz 


signal_l; 


s_b « 


signal_2 ; 


s_c = 


signal 3; 


out_a = 


target l; 


out__c = 


target 2; 


factor «= 


1.5; 


scale «= 


1.0 /4160.0; 



/* noise canceller initialization at time t=0 */ 

nc[0] . berr = o.O; 
nc[0). Gamma = 1.0; 



for(m=0; m<NC_CELLS; m++) { 

nc[m].err_a = o.O; 

nc(m j ,err_b = o.O; 

nc(m].Roh_a ^ o.O; 

nc(m].Roh_c = o.O; 

nc[m). Delta « o.O; 

ncfmj.Fsvsqr = 0.00001; 

nc[m].Bsvsqr « 0.00001; 



} 
} 

for(k=0; k<passes; k++) { 



END INITIAXIZATION 



contraction ~ FALSE; 
for(m=0; m< NC_CELLS; m++) { 

nc(m].berrl « ncfmj.berr; 

nc[m] .Bswsqrl = nc[m J . Bswsqr; 



/* Update delay elements 



noiseref = factor * log(1.0 - (*s_a) * scale) 

- log(1.0 - (*s_b) * scale) ; 

nc{0).err_a « log(1.0 - (*s_b) * scale); 

nc(0].err_b = log(l.o - (*s_c) * scale); 



++s_a ; 
++s_b; 
++s_c; 

nc{0] . f err 
nc(0] .berr 
nc[OJ .Fsvsqr 
nc[0] .Bswsqr 



noise_ref ; 
noise_ref ; 

LAMBDA * nc{0). Fsvsqr + noise_ref * noise_ref; 
nc{0) .Fsvsqr; 



/* Order Update */ 

for(n°l;( n < NC_CELLS) 4t (contraction — FALSE); n++) { 

/* Adaptive Lattice Section */ 

m « n-1; 
ii« n-1; 



nc[n]. Delta *■ 
nc[mj. Delta +■ 
Delta_sqr 

nc(n], f ref 
nc(n) .bref 

nc(n) . f err 
nc(n) .berr 



LAMBDA; 

nc(m].berrl * nc(m].ferr / nc(m). Gamma ; 
= nc(m). Delta * nc [ m) . Delta ; 

~ -nc(m] .Delta / nc[ m] . BsvBqrl ; 
= -nc(m). Delta / nc(m] .Fsvsqr ; 

= nc(m].ferr + nc(n).fref * nc[m).berrl; 
= nc[m].berrl + ncfnj.bref * nc[m).ferr; 



nc[n). Fsvsqr «= ncfmj. Fsvsqr - Delta_sqr / nc(m) .Bsvsqrl; 
nc[n). Bsvsqr = nc(m] . Bsvsqrl - Delta_sqr / nc[m) .Fsvsqr ; 

iff (nc[n] .Fsvsqr + nc[ n] . Bsvsqr) > 0.00001 | | (n < 5) ) { 

nc(n) .Gamma = nc[m). Gamma - nc(ra).berrl * nc(m).berrl / nc[m) . Bsvsqrl ; 
if (nc(n) .Gamma < 0.05) nc(n). Gamma «= 0.05; 
if (nc[n) .Gamma > 1.00) nc[n). Gamma =1.00; 

/* Joint Process Estimation Section */ 

nc[m].Roh_a *= LAMBDA; 

ncfm].Roh_a +«= nc(m].berr * nc(m).err_a / nc[m). Gamma ; 
nc(m).k_a «= nc(m).Roh_a / nc(m] . Bsvsqr ; 
nc(n].err_a *= nc(m].err_a - nc(m).k_a * nc(m].berr; 

nc[m).Roh_c *«= LAMBDA; 

nc[m).Roh_c +«= nc(m).berr * nc[m].err_b / nc[m). Gamma ; 



nc[m) .k_c 
nc[n) . err_b 



nc[n] .Roh_c 
nc[B) . err_b 



/ nc[m] .Bsvsqr; 
- nc{m].K_c * nc[m).berr; 



) 

else { 

contraction = TRUE; 
for(i=n; i<NC_CELLS; i++) { 

0.0 ; 
0.0; 
0.0; 
0.0; 
0.0; 

0. 00001; 
0. 00001; 
o.ooooi; 



nc[i} .err_a 
nc[i) .Roh_a 
nc[i) . err_b 
nc(i) .Roh_c 
nc[i) .Delta 
nc[ i) .Fsvsqr 
nc[i) .Bsvsqr 
nc(i) .Bsvsqrl 



} 

} 

> 

*out_a++ ■= (int)( (-exp(nc(ii) .err_a) +1.0) / scale) 
*out_c++ « (int)( (-exp(nc(ii).errb) +1.0) / scale) 



> 

) 

/* 



Least Square Lattice 



